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1.  INTRODUCTION 


Continuous  gun  tiring  elevates  the  barrel  temperature,  producing  several  adverse  effects 
on  system  performance.  Accuracy,  and  hence  lethality,  Is  diminished  with  repeated  firings 
due,  in  part,  to  thermal  distortion  of  the  barrel.  Thermal  signature,  and  hence  vulnerability  of 
the  firing  platform,  also  increases  with  firing  rats  and  number  of  rounds  fired.  In  addtfon, 
rapid  firing  of  the  gun  increases  the  concern  that  the  chamber  well  temperature  could  oook  off 
a  subsequent  round.  Barrel  wear  also  Increases  with  gun  tube  temperature.  To  investigate 
the  magnitude  of  these  effects  and  seek  ways  cf  mitigating  their  detrimental  effect  on  the 
overall  gun  system,  the  U.S.  Army  has  embarked  on  a  comprehensive  thennal  management 
program.  Computer  modeling  is  an  integral  part  of  thie  effort,  as  witnessed  by  its  use  in  the 
numerous  reports  being  published  on  these  subjects  (Artus  and  Hasenbein  1989;  Bundy,  to 
be  published;  Chandra  and  Fisher  1989a,  1989b;  Rapp  1990;  Talley  1989a.  1989b). 

Numerical  modeling  Is  the  most  common  approach.  Each  model,  however,  is  developed 
with  specific  objectives  in  mind,  which  makes  its  application  unique.  For  example,  in  the 
multiple-round  cook-off  studies  of  Chandra  and  Fisher  (1989a,  1989b),  emphasis  is  placed  on 
accurately  modeling  the  barrel  wall  temperature  in  the  combustion  chamber  over  the  relatively 
short  time  of  the  combustion  and  blowdown  cycle,  with  no  attention  given  to  the  post¬ 
blowdown  external  cooling  effects  between  rounds  (which  are,  admittedly,  small  In  the  gun 
chamber  region).  Similarly,  In  the  single-round  barrel  wear-type  studies  of  Talley  (1989a, 
1989b),  attention  is  focused  on  the  bore  surface  temperature  over  the  first  half-second  after 
firing.  On  the  other  hand,  in  the  multiple-round,  full-barrel,  first-cut  thermal  surveys  of  Artus 
and  Hasenbein  (1989)  and  Rapp  (1990),  less  detailed  (time-averaged)  propellant  heat  Input  is 
used,  with  some  consideration  given  to  external  cooling  (through  the  use  of  free  and  forced  air 
convection  assumptions).  The  multiple-round,  full-barrel  temperature  model  of  Bundy  (to  be 
published)  uses  experimentally  measured  external  cooling  rates,  but  it  is  empirical  In  nature 
and  thus  limited  to  the  range  of  operating  conditions  upon  which  it  was  developed. 

We  seek  to  establish  a  full-barrel  temperature  modeling  capability  for  the  M1A1, 120-mm 
M2 56  tank  gun.  The  method  of  solution  will  be  finite-difference  based  and  similar  to  that  of 
Chandra  and  Fisher  (1989a,  1989b).  External  cooling  rates  will  be  based  on  experimental 
data  obtained  for  this  particular  (field-configured)  gun  system  (with  thermal  shroud,  bore 
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evacuator,  muzzle  reference  system  collimator,  and  standard  M1A1  raool  mount  system). 
Eventually,  we  wish  to  modtiy  the  programming  to  simulate  passive  and  active,  Internal  and 
external  barrel  cooling  effects,  with  the  long  range  goal  of  developing  a  capebfllty  to 
investigate  the  feastoility  of  various  gun  barrel  cooling  devtoes. 

We  intend  to  develop  this  model  Incrementally,  documenting  It  in  a  series  of  reports, 
beginning  herewith.  This  report  concerns  barrel  heating  and  cooling  In  one  dimension  (racial). 
Nevertheless,  it  will  be  possible  to  Illustrate  the  importance  of  aocurats  round-to-round  heat 
input  data  in  successfully  predating  the  temperature-time  history  within  the  barrel. 

The  two  problems  of  determining  the  flow  in  the  be  re  and  the  heating  in  the  barrel  are 
coupled  in  that  both  involve  the  temperature  at  the  inner  wall  of  the  barrel.  In  principle,  an 
iterative  procedure  should  be  applied  between  the  two  problems;  however,  that  is  frequently 
not  practical.  We  shall  perform  only  the  first  approximation  here,  assuming  that  the  flow 
problem  'drives*  the  heating  problem,  in  this  case,  the  interior  flow  equations  usually  contain 
an  approximation  for  convective  heat  loss  to  the  barrel. 

We  shall  be  following  the  lead  of  previous  investigators  In  assuming  that  heat  conduction 
in  the  axial  direction  may  be  neglected  relative  to  that  in  the  racial  direction  (e.g.,  see  Comer 
(1950],  Engineering  Design  Handbook  (1965],  Heiney  (1979]).  Furthermore.  If  we  ignore  the 
effects  of  gravity  on  convective  heating  and  cooling,  then  the  problem  is  axisymmetf  c,  and 
there  is  no  azimuthal  conduction  in  a  transverse  plane.  Thus,  we  begin  by  treating  only 
unsteady  radial  heat  flow  through  the  annular  barrel.  Figure  la  shows  a  longitudinal  cross 
section  of  the  barrel  of  an  M2 56  120-mm  gun  (vertical  scale  expanded).  We  should  expect 
that  our  model,  at  this  stage,  would  be  most  applicable  at  locations  away  from  the  comers. 

To  keep  the  model  tractable,  we  make  some  further  approximations.  We  neglect  any  heat 
addition  to  the  barrel  produced  by  shock  wave  heating  ahead  of  the  projectile  and  by  friction 
heating  caused  by  the  passage  of  the  projectile.  We  also  omit  thermal  expansion  of  the  metal 
barrel. 
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Figure  la.  Lengthwise  Cross  Section  of  M2 56  120-mm  Gun  Barrel  (Stretched  Vertical  Scale). 

Our  approach,  applicable  to  any  gun  barrel,  is  to  calculate  a  numerical  solution  to  the 
boundary-value  problem  that  simulates  the  heat  transfer  processes.  This  calculation  is  carried 
out  by  the  Crank-Nichoison  implicit  finite  difference  method  (e.g.,  see  C)zisik  [1968, 402]). 

We  shall  perform  computations,  in  particular,  for  a  120-mm  M256  tank  gun  and  compare 
results  with  experimental  data.  Additionally,  we  shall  conduct  a  study  relevant  to  the  simplified 
model  of  Bundy  (to  be  published),  which  provides  a  rapid  approximate  determination  of  the 
heating  of  the  120-mm  M256  tank  gun  by  repeated  firing. 

2.  FORMULATION  OF  THE  PROBLEM 

We  state  our  problem  in  terms  of  the  following  cylindrical  coordinates:  r,  0,  and  z*  The 
radial  coordinate,  r,  is  zero  on  the  axis  of  the  gun  tube  (z-axis)  and  varies  from  rt  to  rv  the 
concentric  radii  of  the  inner  and  outer  walls  of  the  barrel,  respectively  (Figure  1b).  As  stated 
in  Section  1 ,  the  azimuthal  angle,  0,  does  not  enter  the  problem.  The  axial  coordinate,  z.  Is 
taken  to  be  zero  at  the  gun’s  breech.  The  barrel  temperature,  T  (r,  z,  t),  where  t  is  time 


*  Definition*  of  symbols  are  given  in  the  List  of  Symbols  section. 
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Figure  1b.  Transverse  Cross  Section  of  Gun  Barrel. 

measured  from  the  initiation  of  the  first  round,  is  determined  by  the  following  differential 
equation  of  heat  conduction  for  a  stationary,  homogeneous,  isotropic  solid  with  no  internal 
heat  generation  (Ozisik  1968): 

pcp  8779/  =  div  [k  grad  T],  (1) 

Here  the  density,  p,  and  the  specific  heat.  of  the  metal  are  taken  to  be  constant;  *  the 
quantity  k  is  the  thermal  conductivity  of  the  metal.  With  the  simplifications  of  Section  1  and 
the  further  assumption  that  k  is  constant,  Equation  (1)  becomes  the  Fourier  Equation,  namely, 

92778r 2  +  (1//)  dT/dr  *  (1/cx) 8779/,  (2) 

where  a  ■  k/(pc„)  is  the  thermal  diffusivity. 

We  now  introduce  the  following  notation:  7^  -  ambient  temperature  of  the  atmosphere, 

7) »  barrel  temperature  at  the  interior  wall  r  •  r, ,  and  T0  -  barrel  temperature  at  r  -  r0.  We 
assume  an  initially  constant  temperature,  7^,  everywhere.  Thus,  the  initial  condition  at  a 
given  z  is 


*  We  note  that  cp  and  k  are  actually  function#  of  temperature;  the  effect  of  temperature  dependence 
will  be  studied  In  a  future  investigation. 
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(3) 


T(r,z,o)  =  t  *  0,  r,£  r<,  rQ. 

The  boundary  conditions  at  the  inner  and  outer  wails  are  obtained  by  equating  the  rate  of 
heat  transfer  from  the  surrounding  gas  at  the  wall  to  the  rate  at  which  heat  flows  from  the  wall 
into  the  barrel.  The  boundary  condition  at  the  inner  wall  is 

k  BT/dr  -  hgT=  -  hg  Tg  r=r„t>0,  (4) 

where  Tg  (t ,  z)  Is  the  cross-sectional  average  temperature  of  the  combustion  products  in  the 
bore  at  time  t  and  location  z,  and  hg  (t ,  z)  is  the  coefficient  of  heat  transfer  between  these 
products  and  the  inner  wall  of  the  barrel  (e.g.,  see  Ozisik  [1968,  8-9]).  In  our  model,  Tg  and 
hg  will  be  assumed  known  for  any  t  and  z  and  thus  constitute  input  to  the  problem  (see 
Section  3). 

The  boundary  condition  at  the  outer  wail  is 

k  BT/dr  +  hooT  =  /to Too  r  =  rQ,  t  >  0,  (5) 

where  the  constant  is  the  coefficient  of  heat  transfer  between  the  barrel  at  r  *  r0  and  the 
surrounding  atmosphere. 

A  computational  difficulty  arises  at  the  start  of  the  ballistic  cycle  due  to  the  local 
temperature  variation  near  the  inner  wall.  This  problem  is  circumvented  by  introducing  a 
transformed  variable  4, 

r-r(%)  (0  5  5  si).  <6> 

so  that  the  constant  increment  A  4  will  bunch  the  nodal  points  closely  together  near  the  inner 
wall  but  spread  them  out  away  from  there. 

We  define  our  transformation  in  the  following  two  steps: 

4  =  y4  +  (1  -  y)  4P  (o<y  <  1,  P  >  2) 

r  =  D  4  +  r,f  (7) 
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where  D  «  r0  -  r, ,  the  barrel  thickness,  and  y  and  0  are  chosen  constants.  We  have  used 
Y  -  0.092,  p  -  2.25.  Note  that  r  -  r, ,  r0  correspond  to  4  -  0,1.  Then  Equation  (2)  transforms 
to 

(8) 


(9) 


and  where  ()'  =  (/(  )ld%\  the  formulas  for  C/  and  are  given  in  Appendix  A. 

Equation  (4)  transforms  to 

k  zm%  -  Dh'  (%)  t  -  oh,  C'  (4)  r„  4  =  o,i>  o,  (io) 


and  Equation  (S)  transforms  to 

*97/94  +  DhJL (4)7=  0/ipo^  (4)7oo  4  =  1,f  >  0.  (H) 

Then  Equations  (8),  (10).  (11),  and  (3)  comprise  the  actual  problem  to  be  solved  numerically. 
In  the  final  printout,  we  revert  to  use  of  the  independent  variable,  r,  given  by  Equation  (7). 

In  the  present  model,  we  assume  that  NR  rounds  are  fired  at  the  constant  interval,  t, , 
between  rounds.  The  input  functions  Tg  and  hg  are  taken  to  be  periodic,  over  period  f, ,  at  all 
stations  on  the  barrel.  The  subsequent  barrel  cooling  after  NR  rounds  can  also  be  calculated. 


BT/dt  -  (a ID2)  (/;  (4)  d2T/ 9  42  +  f2  (4)  97794]  «  G  (4.0. 

where 

/.  (4)  =  - 1 - 

IC'<4))2 

#/« .  pirn  m 

2  “  D  4(4)  ♦  r,  [C/(4)l3  ' 
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3.  INPUT  DATA 


3.1  Bore  Temperature  and  Convective  Heat  Transfer.  As  stated  previously,  the 
temperature,  Tg  (t,  z),  and  the  convective  heat  transfer  coefficient,  hg  (t,  z),  of  the  gas  In  the 
gun  bore  are  provided  as  input  at  the  inner  wall.  The  Tg  and  hg  values  are  obtained  from 
calculations  modeling  the  flow  in-bore  during  a  firing  cycle.  There  are  a  number  of  models  of 
the  interior  ballistics,  with  varying  degrees  of  realism  in  simulation.  Most  of  these  compute 
"average"  values  for  flow  variable  functions  by  ordinary  differential  equations.  The  NOVA 
code  (Gough  1980),  however,  includes  axial  variation  in  the  computation  of  the  flow  variables. 

It  is  the  Veritay  modification  of  this  model  (Chandra  and  Fisher  1989a,  1989b)  that  supplies  Tg 
histories  at  eight  chosen  stations  along  the  barrel.  The  values  of  ha  furnished  by  the  Veritay 
code  are  obtained  from  the  flow  variables  on  the  basis  of  the  method  of  Stratford  and  Beavers 
(1961). 

In  the  simplest  simulation,  one  would  model  the  in-bore  flow  for  a  single  round  out  to 
about  1 ,000  s  and  would  repeat  as  many  times  as  needed  only  that  portion  of  the  history  that 
covers  the  time  interval  0  5 1  s  t,.  In  practice,  there  are  factors  that  complicate  the  flow  model, 
such  as  loading  a  new  round  into  the  chamber.  These  factors  will  be  neglected  here. 

A  difficulty  arises  with  use  of  the  current  NOVA  modei  in  that  the  numerical  process  stops 
operating  at  some  time  during  the  blowdown  after  the  projectile  exit  from  the  muzzle.  At 
present,  the  only  means  of  extending  the  Tg  and  hg  curves  is  extrapolation.  Each 
extrapolation  curve  will  be  required  to  have  the  same  position  and  slope  as  the  input  curve  at 
some  time,  f8,  and  to  approach  asymptotically  the  ambient  value  of  the  variable  as  «. 

We  employ  an  exponential  extrapolation  for  1 2  tg,  where  the  break  point  tg,  is  taken  to  be 
0.065  s.  The  time  derivative,  T,  is  approximated  at  t  -  tB  by  a  difference  quotient 
[Tg  (tB  +  Atg)  -  Tg  (tB  -  Atg)j/2  Atg,  where  A tg  is  the  time  increment  in  the  Tg  tables  and  TB  *  Tg  (tg). 

T,  (t)  -  T.  =  (T,  -  TJ  exp  [  f,  (f  -  l,)l(T,  -T  )].  <12) 

A  similar  formula  is  used  for  hg(t).  Figures  2a  and  2b  show  representative  Tg  and  hg  histories 
at  two  stations  on  a  120-mm  gun  barrel.  It  is  seen  that  Tg  and  hg  remain  constant  until  the 
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Figure  2b.  Convective  Heat  Transfer  Coefficient  at  Inner  Wall 
Locations. 


base  of  the  projectile  passes  the  given  station  at  time  t  -  ^ .  At  this  time,  these  variables  rise 
suddenly,  then  they  decrease  more  gradually,  with  hg  decaying  significantly  faster  than  Tg . 

3 2  Gun  Properties  and  Ambient  Conditions.  All  the  computations  reported  here  were 
performed  for  the  case  of  an  M256  120-mm  tank  gun  firing  a  DM13  round.  The  dimensions  of 
the  gun  barrel  are  shown  in  Table  1 . 


Table  1.  Dimensions  of  M256  Gun  Barrel 


z 

(mm) 

2r, 

(mm) 

z 

(mm) 

2r0 

(mm) 

z 

(mm) 

2r0 

(mm) 

0.0 

159.1 

000.0 

270.0 

2,788.0 

225.0 

61.0 

157.6 

237.0 

270.0 

2,804.0 

225.0 

486.0 

157.6 

238.0 

309.9 

2,805.0 

216.0 

555.0 

120.8 

800.0 

309.9 

2,868.0 

216.0 

805.0 

120.1 

1,000.0 

270.0 

2,898.0 

178.0 

5,300.0 

120.1 

1,731.0 

270.0 

3,415.0 

178.0 

1,732.0 

250.0 

3,445.0 

216.0 

1,863.0 

250.0 

3,520.0 

216.0 

1,864.0 

240.0 

3,521.0 

210.0 

1,900.0 

240.0 

3,560.0 

210.0 

1,901.0 

225.0 

3,561.0 

171.0 

1,940.0 

225.0 

5,030.0 

161.0 

1,950.0 

215.0 

5,031.0 

155.0 

2,762.0 

196.0 

5,300.0 

155.0 

The  values  of  properties  of  the  gun  barrel  metai  are  tak?n  to  be 

Cp  =  469.05  J/(kg  deg  K)* 
k  ^  38.07  J/{m  s  deg  K) 
p  -  7,827.0  kg/m3. 

*  The  above  value  of  Cp  was  measured  by  Joseph  Cox,  Benet  Weapons  Laboratory  (1 990),  for  M256  gun 
barrel  steel  at  295  K;  the  values  for  k  and  p  were  obtained  from  Talley  (1 989b)  for  4335  steel;  the  value  of 
hm  was  obtained  from  experiments  conducted  by  Bundy  on  a  shrouded  M256  barrel. 


The  ambient  condition  constants  are 


T  m  294.8  K 

oo 

hm  -  12.0  J/(m*s  deg  K). 

4.  FINITE-DIFFERENCE  CALCULATION 

In  the  Crank-Nichoison  method  employed  here,  all  derivatives  (except  5  derivatives  at  the 
walls)  are  approximated  with  central  difference  expressions.  At  the  walls.  \  derivatives  are 
obtained  with  three-point  formulas.  As  a  consequence,  at  each  time  step  the  temperature 
profile  is  determined  by  solving  a  set  of  NI+ 1  simultaneous  linear  equations,  where  Nl  is  the 
number  of  intervals  formed  by  the  nodal  (or  grid)  points  between  the  inner  and  outer  walls. 

The  constant  £  increment  is  given  by  A 5  «  1/NI.  Figure  3  shows  a  diagram  of  the  grid  scheme 
for  approximating  Equation  (8)  at  point  P.  The  temperature  is  determined  at  time  in  terms 
of  the  temperature  at  the  previous  time  step,  tt «  /  At  and  the  boundary  condition  at  /  -  tM, 
where  At  >s  the  time  increment. 

We  begin  with  the  boundary  conditions.  At  4  *  0,  Equation  (10)  is  approximated  by 

[k/ (2D  AS,  C'  (0)11  (37T  -  4 7j*’  <-  T,'”)  =  T,  -  h,  T,"'.  (13) 

Similarly,  at  %  •  1,  Equation  (11)  is  approximated  by 

[*/(2D  AS,  C'  (1)))  (Tg,  -  4 7*’  -  3Tnm)  =  ft.  T.  -  ft.  T&\ ,  04) 

where  the  subscript  index  denotes  the  nodal  point,  and  the  superscript  index  denotes  the  time. 

For  the  interior  points,  2  s  /  s  Nl,  the  approximation  of  Equation  (8)  reduces  to 

7^  -  (Atl2)G^  =  T/  +  (At/2)  Gj,  (15) 
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Figure  3.  Grid  Diagram  for  Numerical  Solution. 
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where  the  right  hand  side  is  now  Known.  The  formula  for  Gj  ( and  Gj*x )  contains  3T/3$  and 
3^/3$*.  These  are  approximated  by 

(3 T/Wr(Thi  -  Th)/( 2A4) 

0*773  V),  -  ♦  Th  -  2  ry)/(A4)a.  (16) 

Then  Qi  for  both  #  A t  and  (i+ 1)  At  is  given  by  the  linear  expression 

g,  .  g,  <y  r;„  .fttyij.ft  (V  V  (i7) 

Formulas  for  g„  p„  and  g9  are  given  in  Appendix  A.  Thus,  Equations  (13),  (14),  and  (15) 
provide  us  with  a  set  of  linear  equations  for  the  temperature  at  NI+1  points  at  time 
t-(to1)At: 

52  A,.  1?'  -dr  u  -  1,  2.  ...  Nl  *  1).  (18) 

n  •  1 

The  coefficients  and  dj  are  given  in  Appendix  A.  A  standard  FORTFIAN  routine  is  applied 
to  soive  Equation  (16);  we  have  in  most  cases  used  Nl  -  125. 

There  is  a  problem  in  choosing  At  because  there  are  essentially  two  time  scales—  (1)  the 
duration  of  the  firing  (roughly  100  ms)  and  (2)  t,,  the  time  between  firings  (usually  5  s  or  more 
when  firing  large  guns).  The  At  should  be  sufficiently  small  to  resolve  the  phenomenon  in 
case  1  but  should  be  larger  in  case  2  to  save  time  in  computation.  The  program  contains  a 
subroutine  prescribing  At  as  a  function  of  f  within  a  firing  cycle  (see  Appendix  B). 

The  coefficients  in  the  heat  conduction  equation  (Equation  8)  are  assumed  to  be 
independent  of  t  Thus,  only  a  single  iteration  is  required  to  obtain  the  solution  to  the  finite- 
difference  equations.  The  Crank-Nicolson  method  is  stable  for  all  values  of  At,  and  there  are 
no  restrictions  on  the  relative  sizes  of  At  and  Alj. 
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5.  ACCURACY  CHECKS 


At  every  time  step,  the  following  integral  quantities  are  numerically  evaluated: 


t 


Og  (t)  -  2  n  r,  j  hg  [Tg  (q)  -  T  M)]  dr\ 

0 

(19) 

i 

o.  (/)-2nr0J/i.[r.  -  T </•„!,)]  dr\ 

0 

(20) 

r. 

Qb  (t)  -  2  71  pcp  f  r  [T(r,t)  -  7J  dr. 

(21) 

r, 


Q,  m  (Qg  +  Q'j  is  the  quantity  of  heat  per  unit  length  of  barrel  that  has  entered  the  banal 
through  the  inner  and  outer  walls,  and  Q„  Is  the  increase  in  the  quantity  of  heat  per  unit  length 
within  the  barrel  since  t  -  0,  A  necessary,  but  not  sufficient,  condition  for  accuracy  Is  that 
Qbm  Qt.  With  the  assumption  of  no  errors  In  the  code,  poor  agreement  generally  indicates 
that  Af  and/or  a£  should  be  decreased.  Empirical  studies  of  temperature  output  vs.  At  and  A£ 
were  additionally  used  as  guides  in  choosing  numerical  parameters. 

Figures  2a  and  2b  show  that  Tg  and  hg  experience  sharp  jumps  at  points  along  the  bore 
just  as  the  base  of  the  projectile  passes  those  locations.  This  produces  a  timewise 
discontinuity  in  the  inner  wall  boundary  condition.  To  study  the  'damage*  that  such  a 
singularity  might  cause  in  our  numerical  output  we  solve  a  simpler  problem,  both  by  our 
numerical  scheme  and  by  an  approximate  analytical  method  that  is  not  affected  by  the 
singularity.  The  approximate  method,  which  we  designate  the  “thermal  layer*  method,  is 
described  in  Chapter  7  of  Ozisik  (1968).  It  is  applicable  only  at  very  early  times;  an  outline  is 
provided  in  Appendix  C.  In  this  simpler  problem,  which  differs  from  the  main  problem,  the 
following  conditions  hold: 


7-0  for  f  s  0,  ry  £  r  £  r0 

Tg  -  Tgg,  ha  -  />„  for  />  0,  (22) 


where  and  are  chosen  constants;  in  addition,  Equation  (4)  applies,  and  Tm  -  0. 

Figure  4  shows  the  temperature  histories  at  the  Inner  wait  determined  by  the  two  methods 
for  a  given  set  of  parameters.  The  good  agreement  between  the  two  outputs  is  an  indication 
that  the  finite  difference  method  is  not  seriously  affected  by  impulsive  changes  of  variables  in 
the  bore. 

6.  COMPUTATIONAL  RESULTS 

6.1  Enemy  Considerations.  The  first  topic  of  interest  is  the  energy  transferred  from  the 
bore  to  the  gun  barrel.  Measurements  (Talley  1989b;  8rosseau  et  al.  1982)  have  been  made 
of  Qa,  the  total  heat  transferred  per  unit  area  of  the  inner  wall  for  a  single  round.  The  quantity 
Ga  is  related  to  (^,  computed  here,  by  the  relation  QA  -  <V(2x  r,).  Figure  5  shows  computed 
values  of  QA  at  five  locations  on  the  barrel  of  an  M256  120-mm  gun  that  has  fired  a  DM13 
round.  Also  shown  are  experimental  results*  of  Tally  (1989b)  and  Brosseau  et  al.  (1982). 

The  largest  discrepancy  is  roughly  19%  at  z  -  5.18  m;  reasonably  good  agreement  is  obtained 
for  the  other  three  comparisons. 

Figure  6  provides,  along  with  gas  temperature  and  heat  transfer  coefficients,  a 
representative  history  of  heat  delivery  to  the  barrel  for  a  single  round.  The  location, 
z  -  2.78  m,  is  approximately  halfway  along  the  length  of  the  barrel.  It  Is  seen  that  at  this 
location  about  66%  of  the  heat  Is  transferred  after  the  projectile  has  left  the  gun  at  t  -  7.5  ms. 
Practically  ail  the  energy  has  been  deposited  by  t  -  40  ms  even  though  the  temperature  of  the 
gas  in  the  bore  is  still  decreasing  significantly.  Figure  6  indicates  that  the  reason  for  this  is 
the  decay  of  the  heat  transfer  coefficient  (h,)  to  a  very  small  value  by  this  time. 

*  The  numbers  of  Talley  end  Brosseau  were  adjusted  by  multiplicative  lectors  In  order  to  match  our  use  o! 
c>  -  489.05  joules/(ms  K)  as  the  specific  heat  of  the  metal. 
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Figure  6.  Heat  Gain  of  Gun  Barrel  During  Firing  of  Round  at  z  =  2.78  m. 


62  Speod  of  Heat  Penetration.  It  takes  a  finite  time  after  firing  for  a  detectibie 
temperature  rise  to  be  measured  at  a  given  point  inside  the  barrel  or  on  the  outer  wall.  To 
obtain  an  estimate  of  the  speed  of  heat  penetration,  we  define  the  history  of  a  "heat  pulse"  as 
the  locus  of  the  temperature  T  «  T_  +  0.5  K  in  the  r,  t  plane.  (0.5  K  is  a  change  that  can  be 
measured  by  thermocouples.) 

Figure  7  shews  records  of  the  heat  pulses  at  three  stations  on  the  barrel.  The  three 
curves  coincide  for  the  durations  of  the  z  «  3.67  m  and  z  -  5.18  m  pulses.  The  times  to 
reach  the  outer  wail  are  approximately  19.5  s,  4.3  s,  and  1.8  s  at  z  =  2.78  m,  3.67  m  and 
5.18  m,  respectively.  For  the  most  part,  the  velocity  decreases  with  time.  An  exception 
occurs  at  z  *  2.78  m  (r0  »  0.1 1  m),  where  the  pulse  accelerates  slightly  after  about  14  s.  A 
test  calculation  with  roa0.16matz>  2.78  m  produced  a  similar  behavior,  the  acceleration 
beginning  at  about  72  s. 

6.3  Outer  Wall  Temperature:  Slow  Rate-of-FIre.  Figures  8a  and  8b  show  the 
temperature  histories  for  several  radial  stations  at  two  axial  locations  computed  for  a  slow  rate 
of  fire,  namely,  one  round  every  4  min.  It  is  observed  that  for  each  round  there  is  a  time  after 
firing  when  all  the  curves  practically  coalesce,  implying  the  attainment  of  nearly  constant 
temperature  across  the  barrel.  This  time  has  been  called  the  "soak-out  time."  It  is  also 
roughly  equal  to  the  "rise  time,"  t,,  namely,  the  time  (for  a  given  round)  when  T0  attains  its 
maximum.  We  shall  henceforth  just  use  the  "rise  time"  expression.  One  cannot  designate 
precisely  a  value  for  t,,  but  one  can  estimate  an  approximate  value  from  the  curves.  See, 
e  g..  Figure  8a,  first  round.  The  rise  time  value  seems  to  stay  constant  from  round  to  round  at 
a  given  station,  but  it  changes  with  barrel  thickness.  Titus  t  -  100s  at  z  -  2.78  m  and  -  14  s 
at  z  »  5.18  m. 

The  temperature  behavior  described  above  has  been  found  to  be  fairly  general. 
Henceforth,  "slow"  and  "fast"  rates-of-fire  will  refer  to  situations  in  which  t,  <  t,  and  t,  <  t,, 
respectively,  where  t,  is  the  interval  between  rounds. 

Figure  9  presents  a  comparison  between  theoretical  and  experimentally  determined 
(Bundy,  to  be  published)  values  of  rise  of  outer  wall  maximum  temperature  between 
successive  rounds,  3  min  apart,  at  z  =  4.30  m.  The  two  sets  of  points  agree  to  within 
about  10%. 
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..  Slow  Rate-of-Flre  (t| »  4  min):  Temperature  Histories  at  z  «  2.78  m. 


Figure  8b.  Slow  Rate-of-Fire  (t  «  4  min):  Temperature  Histories  at  z  =  5.18  m. 
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6.4  Outer  Wall  Temperature:  Fast  Rate-of-Fire.  Figure  10  shows  the  Inner  and  outer  wall 
temperature  variations  computed  at  z  -  4.3  m  (1  m  from  muzzle)  for  a  15-round  burst  of 
three  rounds  per  minute.  The  increase  In  T0  between  the  beginnings  of  the  fifth  and  fifteenth 
rounds  is  77.4  K.  This  number  compares  favorably  with  the  corresponding  experimental 
figure,  approximately  75.7  K.  read  from  the  graph  in  Figure  13  of  Bundy  (to  be  published). 

7.  RELATION  TO  ANALYTICAL  -  EMPIRICAL  MODEL 

We  examine  some  aspects  of  our  simulation  in  relation  to  the  simplified  analytical  model  of 
Bundy  (to  be  published)  for  gun  tube  heating.  Mathematically,  this  model  requires  only  the 
evaluation  of  formulas.  It  is  also  empirical,  in  that  certain  physical  parameters  employed  are 
determined  experimentally.  Two  advantages  of  this  model  are  its  low  computational  time  on 
the  computer  and  its  ability  to  simulate  a  large  variety  of  firing  sequences.  So  far,  only  those 
parameter  values  that  are  applicable  to  an  M256  120-mm  gun  firing  DM13  rounds  have  been 
used.  Figure  1 1  shows  a  typical  example  of  results  produced  by  this  analysis. 

For  slow  rate-of-fire.  the  model  assumes  that  the  outer  wail  temperature  has  a  sawtooth¬ 
like  variation  with  time,  as  in  Figure  11.  For  each  firing  cycle,  the  left  segment  (always  with 
positive  slope)  extends  from  the  instant  of  fire  through  tp  the  rise  time  interval.  The  right 
segment,  with  uniform  temperature  across  the  barrel,  extends  to  the  next  instant  of  Are.  The 
rise  times  have  been  determined  experimentally  by  temperature  measurements  on  the  outer 
wall  of  the  barrel.  It  can  be  seen  that  this  modeling  is  an  idealization  of  the  heating  pattern 
depicted  in  Figure  8.  In  fact,  rise  times  estimated  from  computations  may  be  compared  with 
measured  values  for  checking  purposes. 

For  fast  rate-of-fire.  the  same  sawtooth  pattern  cannot  be  used  to  simulate  temperature 
time-variation.  Figure  10,  where  t,  / 1,  -  0.72,  demonstrates  this  situation.  Modeling  will 
necessarily  be  more  complex  here.  Details  of  the  treat nent  of  this  case  in  the  Bundy  model 
are  deferred  to  Bundy  (to  be  published).  One  assumption,  however,  is  that  the  heat  input  to 
the  gun  from  each  round  will  be  complete  before  the  next  round  is  fired.  This  assumption  is 
certainly  valid  for  the  conditions  in  Figure  6,  which  is  a  representative  example  of  heat  transfer 
history.  The  heat  input  time  is  about  40  ms,  while  times  between  rounds  are  at  least  of  the 
order  of  several  seconds. 
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Figure  11.  Outer  Wail  Temperature;  Bundy  Model  and  Experimental  Results  (Sequence  of 
4  Slow.  14  Rapid.  13  Stow  Roundel. 

8.  DISCUSSION 


A  limited  number  of  comparisons  for  a  single  gun  system  were  made  in  Section  6  between 
i_  juts  of  the  present  computer  model  and  experimental  data.  The  resulting  favorable 
acumen  Vindicates  that  this  simulation  can  yield  reasonable  predictions  of  gun  tube  heating. 

It  our  intention  to  remove  some  of  the  limitations  of  the  present  model.  Extension  to  the 
two  lensionai  (radial  and  axial)  unsteady  heat  conduction  problem  will  undoubtedly  be 
laborious.  However,  two  less  difficult,  but  nonetheless  significant,  refinements  can  be  made  to 
the  one-dimensional  model.  The  first  is  the  addition  of  a  thin  layer  (0.15  mm)  of  chrome  to 
the  inside  wall  of  the  barrel.  (This  layer  is  platted  onto  the  gun  to  decrease  erosion.)  The 
second  modification  is  to  introduce  more  accurate  input  for  the  physical  parameters  in  the 
conduction  equation  and  boundary  conditions.  Thus,  for  example,  Cp  and  K  are  actually 
functions  of  temperature.  This  refinement  would  render  the  problem  analytically  non-linear. 

It  is  expected  that  future  interior  ballistic  input  will  be  produced  by  an  updated  version 
(XKTC)  of  the  NOVA  code.  This  new  version  will  be  more  "robust*  than  the  present  one,  and 
will  carry  the  blowdown  calculation  further  out  in  time. 
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1 .  From  Equation  (7),  4  -  Y  £  ♦  0  -  Y)  4P-  Then 

<*C /<**■;'  -Y*P(1  -Y)^1 

(f  Urf*  ■  4"  -  p  (p  -  1)  (1  -  Y)  £M  <A>1) 


2.  The  functions  g,(4),  ga(4),  and  g,(4)  occurring  in  Equation  (17)  are  evaluated  by  the 
following  sequence: 

A  «  1/(C02 

4  -  -  a)  -  cm 

01  -  i*Da)[y»A  4)a  -  V(2  A  4)] 

ft*  -(ZaiWA*  4)2 

ft  *  (a/o»)  [//(a  V’  ♦  4/(2  a  4)]. 


(A.2) 


3.  The  coefficients  A ^  and  dt  In  Equation  (18)  are  now  given: 
For  J  -  1, 


For  J  m  NI+1, 


Ai  -  3  *  2  A  4  DC7  (4  »  0)  V* 

A,2  *  -4,  A,3  ■  1. 

ijvf.  i  *  ^(2  A  4),  t  Mj  *  -2/A  4 
=  3/(2  A£)  */>.  OC(4»iy* 


(A.3) 


(A.4) 
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Fo rSiJiNI, 


*A,_,  - -(A*2)Mty  ^  -  D  A  ^ 

-  1  -  (A*2)  Ai  (ty 

*  -  (AW)  Mty-  (A.5) 


All  other  coefficients  \  are  equal  to  zero. 

d,  -  2  A  5  o  C'  «  -  0)  h,  TJk 
dm.,  -  K  T.D?  (5  .  iy*.  (A-#) 


Forays  A// 


d,-  Tf  +  (A W)G/,  where 

<v  !?..♦*  <y»7*  mv7?-" 
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APPENDIX  B: 
TIMESCALE  SUBROUTINE 
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Here,  time  =  f'  will  refer  to  time  within  one  firing  cycle;  f'  =  0  at  the  b  ginning  of  the 
cycle.  Six  constants  are  given:  td,  t,,  tf,  A and  A t£.  Here  4  is  the  uelay  time  for  the 
rapid  rise  in  Ta  and  hg  from  initial  conditions,  and  f,  is  the  time  between  successive  firings. 
The  time  increment  At  (tf)  is  given  by  the  following  function: 


where 


A  t  =  td  0  £  t'  <  ta 

At  =  At{  td  £  t‘  <  t[ 

Af  =  C,  +  C2t'  ti  £  t'  <  £ 

At  =  Af2  fa  ^ 


C2  =  (Afa  -  At[)!(t2  -  t[)  and  C,  =  A t2  -  Cz  t2 


(If  t‘  +  Afa  >  f,,  set  Af  =  t,  -  r). 


A  typical  set  of  values  of  the  parameters  would  be  the  following: 

f,'  =  0.018  S,  fa  =  10.0  S,  Af(  =  0.00025  S,  A f2  =  6.0  S. 
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APPENDIX  C: 

OUTLINE  OF  THERMAL  LAYER  METHOD 
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This  method  (dzisik  1968)  assumes  that  at  early  time  all  the  heat  transferred  from  the 
bore  to  th6  barrel  lies  In  a  thin  layer,  of  thickness  6 (t),  adjacent  of  the  inner  wall.  For  practical 
purposes,  we  can  apply  the  conditions 

T  m  0,  3773 r  ■  0  at  r  ■  r,  *  3(f).  (C.1) 

Assume  the  following  approximate  form  for  temperature  variation: 

T  -  b(t)  [r  -  ( r,  +  5)ja  In  [r/(rt  >  5)].  (C.2) 


This  form  automatically  satisfies  the  conditions  of  Equation  (C.1).  By  applying  Equation  (4)  to 
Equation  (C.2),  we  obtain 

ft  -  ft,c  V{(2*5  +  h9c  *)  In  Mr,  ♦  5)]  -  k  6*/r,}.  <C  3> 

where  and  7^  are  chosen  constants  in  Equation  (22). 


Finally,  we  apply  the  heat-balance  relation  obtained  by  multiplying  both  sides  of 
Equation  (2)  by  r,  integrating  from  r  =  rltor=rl  +  5 (t),  and  applying  Equation  (C.1): 


-  [r  3773r]r  -  (1/a)  (tfcfl) 


(C.4) 


This  leads  to  a  differentia)  equation  of  the  form  db/ctt  -  f(5).  After  some  labor,  a  solution  is 
obtained  in  /  *  f  (&)  form  given  by  the  following  sequence  of  formulas: 

e  -  6 !r,  (C.5a) 


o  -  kl(r,  h^) 


(C.5b) 

* 
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t  -  a  tlrf 


(C.5c) 


c,  ■  3  -  60 


(C.5d) 


»  9o 


(C.5e) 


IcSln 
2  «f 


(i- 


23  ^ 
20 


(1-in 

l2 


(C.5f) 


B  (e)  •  C,  e  ♦  Cj 


(C.5g) 


In  the  above  derivation,  several  approximations  were  made  via  series  expansions  on  the 
assumption  that  e  «  1 . 


t 
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LIST  OF  SYMBOLS 


coefficient  in  linear  equations  for  barrel  temperature,  Equation  (18) 
cp  specific  heat  of  gun  barrel  fioules/kg  K)] 

D  *  r0  -  r,  [m,  mm],  thickness  of  gun  barrel 

dj  coefficient  in  linear  equations  for  barrel  temperature,  right  hand  side  of 
Equation  (18) 

fv  f2  given  functions  of  Equation  (3) 

Gfc.t)  function  for  £  and  /  defined  in  Equation  (8) 

gvg2,  g$  functions  of  f;  in  Equation  (17),  defined  in  Appendix  A 

hg  heat  transfer  coefficient  -  bore  gas  to  gun  barrel  Qoules/fm^  K)] 

/>„  heat  transfer  coefficient  -  gun  barrel  to  ambient  air  Qoules/fn^s  K)] 

j  index  indicating  radial  location  of  a  nodal  point  in  finite  difference  calculation 

k  thermal  conductivity  of  gun  barrel  Oou!es/(m  s  K] 

I  index  indicating  time  at  which  temperature  is  calculated 

Nl  number  of  interv? Js  in  r,  s  r  s  r0  formed  by  the  nodal  points 

NR  number  of  rounds  fired 

Oa  •Q„/(2Krl) 

Qb  increase  in  quantity  of  heat  in  gun  barrel  since  f  -  0,  per  unit  length  of  barrel 
Qoules/m] 

Qg.  Qx  quantities  of  heat  per  unit  length  of  barrel  that  have  entered  the  gun  barrel  through 
inner  and  outer  walls,  respectively,  since  t  *  0  Qoules/m] 

Ot  ■  Qg  + 

r  radial  coordinate  in  transverse  plane  [m,  mm]  (r  ■  0  at  axis  of  gun  bore) 

rA  r0  radial  coordinate  of  inner  and  outer  wails,  respectively,  of  gun  barrel  [m,  mm] 

T  temperature  in  the  barrel  [K] 
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Tg.  Tm  temperatures  in  the  bore  and  ambient  air,  respectively,  [K] 

Tf,  T0  temperatures  at  inner  and  outer  wails,  respectively,  of  gun  barrel  [K] 

Tj  •T£.[J-11ty 

t  time  from  initiation  of  first  round  [s,  ms,  mlri] 

ta  matching  time  for  Tg  and  hg  extrapolation  (Section  3)  [s] 

ta  delay  time  at  given  z  for  rapid  rise  in  Tg  and  hg  [s,  ms] 

tf  time  between  two  successive  firings  [s,  ms] 

tr  rise  time  [s,  ms] 

t'  time  measured  within  a  firing  cycle  [s,  ms] 

two  prescribed  time  values  in  Timescale  subroutine,  Appendix  B  is] 
y  m  r-r,  •  penetration  depth  in  barrel  [m,  mm] 
z  axial  coordinate  (z  -  0  at  breech)  [m] 

a  n  k(p  cj  thermal  diffusivity  of  gun  barrel  ]iyv7s]  Equation  (2) 
p,  y  prescribed  constants  in  transformation  formula,  Equation  (7) 

5  (t)  thickness  of  thermal  layer  In  approximate  method,  Appendix  C 

Ar  radial  distance  between  two  adjacent  nodal  points  [m] 

At  (t)  time  increment  for  calculation  of  temp  rrature  profile  (Equation  15)  [s] 


Atu&tz  two  prescribed  time  Increments  in  the  Timescale  subroutine.  Appendix  B  [s] 
a4  constant  increment  in  %  in  range  0  £  £  s  1 ;  1 1NI 

t,  transformation  variable,  given  in  Equation  (7) 

6  azimuthal  coordinate  in  transverse  plane 

%  transformed  variable,  Equation  (7);  Independent  variable  in  transformed  Fourier 

equation,  Equation  (8) 

p  density  of  gun  barrel  metal  (kg/m3] 
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